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Abstract: We study quantum tunneling for the de Sitter radiation in the planar coor- 
dinates and global coordinates, which are nonstationary coordinates and describe the ex- 
panding geometry. Using the phase-integral approximation for the Hamilton-Jacobi action 
in the complex plane of time, we obtain the particle-production rate in both coordinates 
and derive the additional sinusoidal factor depending on the dimensionality of spacetime 
and the quantum number for spherical harmonics in the global coordinates. This approach 
resolves the factor of two problem in the tunneling method. 
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1. Introduction 

Black hole physics has been a key issue in gravity since Hawking's discovery of black hole 
radiation . Black holes radiate thermal spectra with the temperature determined by the 
surface gravity at the event horizon. One method to get the Hawking radiation is the Bo- 
goliubov transformation between the in-falling state and the out-going state scattered by 
the event horizon [ffl. Recently Parikh and Wilczek have suggested an intuitive interpreta- 
tion of the Hawking radiation as quantum tunneling through the event horizon of particles 
produced from vacuum fluctuations near the horizon [Q] . The tunneling interpretation has 
revitalized the study of Hawking radiation in various black holes. Other phenomena for 
particle production are the Schwinger mechanism || and the expanding geometry [0] . The 
Schwinger mechanism can be understood as quantum tunneling process of virtual pairs 
from the Dirac sea due to the electric field. 

The tunneling method has been mostly applied to static or stationary black holes, 
where the event horizon provides a causal boundary for tunneling of particles from the 
interior to the exterior. Then a question may be raised whether the tunneling process is 
unique to the static or stationary black holes. The de Sitter spacetime, describing the 
geometry of an expanding spacetime, has a few coordinate systems, one of which is the 
static coordinates that exhibit the cosmological event horizon with the Gibbons-Hawking 
temperature ||. The tunneling method has been applied to the static coordinates of de 
Sitter radiation || [7], ||, ^, 10 1. In the static coordinates of the de Sitter spacetime in four 



dimensions [11], |l2| 

r 2 , , n dr 



ds z = -(1 - — )dt z + — +r z dn%, (1.1) 

1 H 



the cosmological event horizon is located at r^s = rjj and has the Gibbons-Hawking 
temperature T^s = l/(27rr#). Volovik used the fluid metric for de Sitter spacetime, which 
is still a stationary coordinates system, and derived the de Sitter radiation via quantum 
tunneling [13]. 



In this paper we study the particle emission via tunneling in nonstationary coordi- 
nates for de Sitter spacetimes. Since Parker has shown that expanding spacetimes also 
create particles 0], the particle-production rate in de Sitter spacetimes has been known 
for many years [ 14 , [l5| . The WKB approximation in complex time has been employed for 



the de Sitter radiation [16]. In the Schwinger mechanism charged particle pairs can also 
be produced by time-dependent electric fields or a constant electric field in time-dependent 
gauge. Hence it is worthy to investigate the de Sitter radiation via tunneling in the planar 
coordinates and the global coordinates, these being nonstationary coordinates. Other mo- 
tivation is the derivation of the Friedmann equation for de Sitter spacetime from the first 



law of black hole thermodynamics 17 ] . 

In expanding spacetimes or time-dependent electric fields, the ingoing positive fre- 
quency solution of a quantum field splits into the outgoing positive and negative frequency 
solutions. The ratio of the outgoing negative solution to the ingoing positive solution de- 
termines the mean number of produced particles [|j . The mean number is approximately 
given by the imaginary part of the Hamilton-Jacobi action, which in turn is determined 
by a contour integral in the complex plane of time |l£|. This approach based on the 
phase-integral approximation ]19| generalizes the tunneling idea to time-dependent electric 
fields [18]. Other field theoretical methods are also attempted in time-dependent contexts 



1 23, pi 



We apply the phase-integral approximation to the planar and global coordinates of 
de Sitter spacetime and obtain the de Sitter radiation. The phase-integral approximation 
for quantum tunneling recovers not only the Boltzmann factor for the particle-production 



rate but also from the Stokes phenomenon 23] the sinusoidal factor that depends on the 



dimensionality of spacetime and the quantum number for spherical harmonics. One contro- 
versial issue in the tunneling method is the factor of two problem, according to which some 
coordinates give the temperature for radiation twice as big as the Hawking temperature 



|24, [25|, 27 1 . We show that the tunneling method in nonstationary coordinates properly 
yields the Hawking temperature and thus resolves the factor of two problem. 

The organization of this paper is as follows. In section 2 we briefly review the Schwinger 
mechanism in the time-dependent gauge and formulate the pair-production rate via quan- 
tum tunneling in the complex time. In section 3 we derive the de Sitter radiation via quan- 
tum tunneling in the planar and global coordinates and explain the origin of sinusoidal 
factor that depends on dimensionality of spacetime and quantum number for spherical 
harmonics. Finally, we summarize the tunneling method in nonstationary spacetimes and 
discuss the problem of the factor of two in section 4. 

2. Schwinger Mechanism in Time-Dependent Gauge 

The quantum electrodynamics (QED) analog for static or stationary black holes is the 
Coulomb gauge, in which the field equation becomes the tunneling problem. The analogy 
between the Schwinger mechanism and the Hawking radiation has been studied in Refs. 



gg , 29, 3C]. In QED the pair-production rate can be found in various methods, such as the 



phase- integral approximation [31, 32, [l8|] and the worldline instanton method [3^, 34]. The 



constant electric field has not only the Coulomb gauge but also the time-dependent gauge, 
in which the field equation becomes either tunneling under the barrier or transmission over 
the barrier. 

We briefly review the phase-integral approximation for pair production in the time- 
dependent gauge []l8|], which is the QED analog for nonstationary spacetimes. The com- 
ponent of a charged scalar field in the gauge potential An (t) satisfies the equation 

0k(t) + <9k(t)0k(*) = O, (2.1) 



where 

Qk(t) = m 2 + ki + (fcy + qA\\(t)f. (2.2) 

Here m and q are the mass and charge of the particle, and and fen are the momentum 
component transverse and parallel to the electric field, respectively. The magnitude square 
of the ratio of the outgoing negative frequency solution to the ingoing positive solution 
is the mean number of produced pairs. In Ref. [18|, the WKB instanton action for the 
Hamilton-Jacobi equation 

2ImS k = i I y/Q\St)dt, (2.3) 

where the contour integral is taken outside of a loop enclosing two roots of Qk(t) in the 
complex plane of time, approximately gives the mean number of produced pairs 

n k = e - 2Im5 *. (2.4) 



The action 2ImS k takes the probability into account, which is equivalent to the detailed 



balance of the emission rate to the absorption rate [28, 21]. 

The instability of the vacuum due to pair production is characterized by the vacuum 
persistence, which is related to the mean number of produced pairs. The vacuum persis- 
tence for bosons is the probability for the in- vacuum to remain in the out- vacuum 35, 36] 



|(out|in)| 2 



-VT£ k ln(l+n k ) 



(2.5) 



where V is the volume and T is the duration. For small pair production % < 1, the decay 
rate per unit volume and per unit time is approximately given by the total mean number 
of produced pairs: 

^ = l-|(out|in)| 2 «^n k . (2.6) 

k 

In the phase-integral approximation for tunneling, it is understood that the emission rate 
is determined by Eq. ( |2.4[ ). 



3. de Sitter Spacetimes 

The de Sitter spacetime has a few different coordinate systems. The (d+l)-dimensional de 
Sitter spacetime has the planar coordinates, 1 

ds 2 = -dt 2 + e 2Ht dx 2 d , (3.1) 

where the units c = H = 1 are used, and also has the global coordinates 

ds 2 = -dt 2 + cosh 2 (Ht)dn 2 d . (3.2) 

The field equation for a scalar field contains the d-dimensional Laplace operator and har- 
monics p7[ 

V 2 d u K (x) = -k 2 u k (x), (3.3) 

where for the planar coordinates V 2 , acts on the d-dimensional Euclidean space and has 
k 2 = k 2 while for the global coordinates acts on and has k 2 = 1(1 + d — 1), (/ = 
0, 1, • • •). Then the harmonic-components of the scalar field, 3>(i, x) = aT d l 2 ^ K it K (x)</> K , 
satisfy 

Mt) + Q K (t)</> K (t) = °> ( 3 - 4 ) 

where 

,. 9 k 2 d.d . /a\ 2 da . 
Q K (t) =m 2 + ---(-- 3.5 

Here a(i) = e Ht for the planar coordinates and a(t) = cosh(Ht) / H for global coordinates. 
We now apply the phase-integral approximation and find the solution of the form 

k(t) = e- lS ^\ (3.6) 

in terms of the Hamilton-Jacobi action 



S K (t) = / VQjfidt. (3.7) 



The imaginary part of S K is responsible for the decay rate and gives the Boltzmann factor 
for the de Sitter radiation, 



r K = |^| 2 = e- 2Im5 *. (3.8) 



First, in the planar coordinates, Eq. (|3.7|) takes the form 



S k (t) = J yV + k 2 e-* m dt, 7 = ym 2 -^L (3.9) 



1 The planar coordinates are also called the oblique coordinates and the global coordinates are the 
horizontal coordinates 



Changing the variable r = e Ht and integrating from a turning point r_ = —ik/j to another 
r + = ik/'y, we obtain a pure imaginary action 



1 dr 
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The analogy with the Schwinger mechanism suggests another way to find the instanton 
action. The production rate of particles may be given by a contour integral III 



2Im5 k = i j V7 2 + k 2 e- 2Ht dt, (3.11) 



in the complex plane of time t. Expanding for large r and choosing a contour exterior to 
two branch cuts r±, the contour integral 



2l ™ s * = i lff + WfT- (3 - 12) 

with the negative sign for residues, leads to the instanton action 

2Im5 k = 2vr-J:. (3.13) 

Second, in the global coordinates, the instanton action of the Hamilton-Jacobi equation 
is given by 



k 2 dr 



There are four turning points, 

e m (a)± = -i(\H/<y) ± i^(XH/^) 2 + 1, 

e m m± = +i(\H/j) ± V( A #/7) 2 + 1- (3-15) 

These are grouped into complex conjugate pairs {t( )+> *(&)-} an d {*(o)-> *(&)+}■ Changing 
the variable to a conformal time, sinh(i?t) = tan(r), the Hamilton-Jacobi action takes the 
form 



Sl(t( a ),t {b) ) = j- 1/1+ 



# V 7 
The integral can be done [38] 



(« 



a) \ 1 + ^ co8(r) 



(3.16) 



Sl(t(a) ,*(&)) = in— +7TA. (3.17) 

Each complex conjugate pair contributes 2Im<Sz = 2iry/H and the crossing integral, for 
instance, from £( )_ to yields the real part ReS; = 7rA. Therefore, the Stokes phe- 
nomenon, Eq. (7) of Ref. [23] (see also Ref. [p^j] ), gives the particle-production rate 

fl, ~ g— 2I mSj(t( a ) + ,t(6)_) _j_ e -2Im5;(i( a )_, t((,) + ) 

+2cos(2Re5 / (i (a) _,i (b) _))e- ImSi ( t w+'*( i ')-)- Im5; ( i W-'*( b )+). (3.18) 



In the asymptotic limit of large action (Si ^> 1 and I 3> 1), we approximate A ~ l + (d— 1)/2 
and obtain 



d. 



^«4 S in 2 (7r(/ + ^)j e -^*, (3.19) 

which is the leading term of particle-production rate in the (d+l)-dimensional de Sitter 
spacetime 



Hl = rinh'fg) ■ (3 - 20) 



The sinusoidal factor is a consequence of the substructure of the Stokes phenomenon and 



explains the absence of particle production in odd dimensions |T5], |39|. In the second 
approach, changing the variable u = cosh(Ht)/H and expanding for large u, the contour 
integral leads to 



2ImS, =ai + = 2*2-. (3.21) 

In summary, we applied the phase-integral approximation to the planar coordinates 
and the global coordinates of de Sitter spacetime and obtained the Boltzmann factor for 
the mean number of produced particles 

n K = e" 2 ^. (3.22) 

The de Sitter radiation has the Gibbons-Hawking temperature T^s = H/2ir. In the global 
coordinates, in addition to the Boltzmann factor, there is a sinusoidal factor that depends 
on the spacetime dimensions and the quantum number for spherical harmonics. That 
particles are not produced in odd dimensional de Sitter spacetimes is a consequence of the 
Stokes phenomenon in the global coordinates. 



4. Conclusion 

In this paper we have extended the tunneling idea for particle production to nonstationary 
coordinates of de Sitter spacetimes. The tunneling method, which provides a physical in- 
tuition to the Hawking radiation of black holes, has been also applied to static coordinates 
of de Sitter spacetimes. It is thus interesting to investigate the tunneling method in non- 
stationary spacetimes. For this purpose we have used the analogy of particle production 
by expanding spacetimes with the Schwinger mechanism by electric fields. 

In QED a constant electric field has the Coulomb gauge or the time-dependent gauge. 
The Coulomb gauge corresponds to a static black hole, while the time-dependent gauge 
corresponds to a nonstationary spacetime. Further, the pair-production rate for a given 
momentum k 



n k = e- 2ImS * 



(4.1) 



is given by the Hamilton- Jacobi action in the complex plane of space or time [18] 



2ImS k = Tij VQUz)dz. (4.2) 

Here the upper sign is for the Coulomb gauge and Q^(z) is the kinematic momentum in 
the electric field while the lower sign is for the time-dependent gauge and Qk(z) is the 



kinematic energy. In this sense Eqs. (4.1) and ( t4.2|) provide a unified tunneling method 



for the Coulomb gauge and the time-dependent gauge. 

We have shown that the tunneling method formulated by Eqs. ( [4,l|) and (L2) also 
applies to the planar and global coordinates of de Sitter spacetimes. The tunneling method 
recovers the Boltzmann factor for de Sitter radiation as shown in section 3 and, to our sur- 
prise, it yields the sinusoidal factor depending on the dimensionality of de Sitter spacetime 
and the quantum number for spherical harmonics in the global coordinates. The absence 
of particle production in the global coordinates of odd dimensional de Sitter spacetimes is 
a consequence of the Stokes phenomenon. 



One controversial issue in the tunneling method is the factor of two problem [24, 25 



|26| , 27\, where the temperature is twice as big as the Hawking temperature. To resolve 



this problem, either the isotropic coordinates [40, 41] or Rindler coordinates p(l have been 



used. Indeed, in these coordinates the tunneling method with the upper sign of Eqs. ( |4.1| ) 
and fl4.2|) yields correctly the de Sitter radiation. However, tunneling method does not 
have such an ambiguity in the planar and global coordinates as shown in section 3. 
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